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Using systematic evolutionary structure searching we propose a new carbon allotrope, pha- 
graphene[fae’graefi:n], standing for penta-hexa-hepta-graphene, because the structure is composed 
of 5-6-7 carbon rings. This two-dimensional (2D) carbon structure is lower in energy than most 
of the predicted 2D carbon allotropes due to its sp^-hybridization and density of atomic packing 
comparable to graphene. More interestingly, the electronic structure of phagraphene has distorted 
Dirac cones. The direction-dependent cones are further proved to be robust against external strain 
with tunable Fermi velocities. 

PACS numbers: 61.46.-w, 73.22.-f, 81.05.U- 

Graphene, as the most stable two dimensional (2D) 
form of carbon, exhibits a number of unusual electronic 
and spintronic properties [1], such as high carrier mobility 
[2] and quantum Hall effects [3, 4]. Its honeycomb lattice, 
which consists of two equivalent carbon sublattices with 
hexagonal symmetry, plays a crucial role in the forma¬ 
tion of the Dirac cones with linear dispersion. Successful 
preparation of graphene in 2004 [5] has inspired further 
searches for other 2D Dirac materials. Among the pre¬ 
dicted 2D Dirac materials [6], are silicene [7], germanene 
[7], silicon germanide monolayer [8], graphynes [9, 10], 

Pmmn boron [11], m-TiB 2 [12], SO-M 0 S 2 [13], etc., while 
only Dirac cones in graphene have been actually con¬ 
firmed experimentally [14]. The stable and robust 2D 
carbon backbones motivate the searches for other Dirac 
carbon allotropes. 

Graphynes, composed of sp- and sp^-hybridized atoms, 
were first proposed in 1987 [15] with various structures 
and high thermal stability. These graphynes can be met- 
alic, semimetalic or semiconducting [16, 17]. Up to 2012, 

Malko et al. [9] reported that a-, /3- and 6,6,12-graphyne 
could have Dirac cones. Huang et al. have tried to derive 
a criterion to explain the physical origin of the existence 
or absence of Dirac cones in graphynes, in which they 
successfully predicted two other graphynes with Dirac 
cones, 14,14,14-graphyne and 14,14,18-graphyne [18]. Af¬ 
ter that, i5-graphyne was proposed as another Dirac car¬ 
bon allotrope, which was demonstrated to be a carbon 
topological insulator superior to graphene [10]. For all 
these graphynes with Dirac cones, a-, /3- and d-graphynes 
exhibit hexagonal symmetry, while 6,6,12-, 14,14,14-, and 
14,14,18-graphynes have rectangular lattices, which sug¬ 


gests that hexagonal symmetry is not a necessary pre¬ 
condition for the presence of Dirac cones. 

However, the sp-hybridization of carbon atoms in¬ 
volved in these graphynes is energetically disadvanta¬ 
geous compared with that of sp^- and sp^-carbon atoms. 
Several carbon allotropes composed of sp^- and sp^- 
carbon atoms have also been predicted to have Dirac 
cones. Liu et al. proposed “T-graphene”, named after a 
buckled carbon sheet with tetrarings, as a 2D Dirac ma¬ 
terial [19]. Although it was demonstrated to be metallic 
in a subsequent work [20] , this opens a possibility to find 
novel 2D Dirac graphene allotropes with multi-member 
rings. In 2014, Xu et al. designed three rectangular al¬ 
lotropes by re-constructing graphene, named as S-, D-, 
and E-graphene. These sp'^- and sp^-carbon networks 
were predicted to be stable and have Dirac cones [21]. 
They also demonstrated that the Dirac cone and the car¬ 
rier linear dispersion are common features of these 2D 
carbon allotropes. 

In this letter, we predict a new carbon allotrope, 
phagraphene, which is composed of 5-6-7 carbon rings. 
This planar carbon allotrope is energetically compara¬ 
ble to graphene and more favorable than other carbon 
allotropes proposed in previous works, due to its sp^- 
hybridization and dense atomic packing. Both density 
functional theory (DFT) and tight binding (TB) model 
confirm the distorted Dirac cone in the first Brillioun 
zone (BZ) for this 2D carbon structure. The direction- 
dependent Dirac cone is further proved to be robust 
against external strain with tunable Fermi velocities. 

Systematic 2D structure searches were performed us¬ 
ing the ab initio evolutionary algorithm USPEX [11, 22- 
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24] with 6, 8, 10, 12, 14, 16, 18, 20, 22 and 24 car¬ 
bon atoms per unit cell. The initial thickness was set 
to zero, since sp- and sp^-carbon allotropes prefer pla¬ 
nar structures at ambient pressure. The newly produced 
structures were all relaxed, and the energies were used 
for selecting structures as parents for the new generation 
of structures. Structure relaxations used the projector- 
augmented-wave (PAW) method [25], as implemented in 
the Vienna ab initio simulation package (VASP) [26, 27]. 
The exchange-correlation energy was treated within the 
generalized gradient approximation (GGA), using the 
functional of Perdew, Burke, and Ernzerhof (PBE) [28] . 

For calculations of specific 2D carbon allotropes in 
VASP, a kinetic energy cutoff of 600 eV was adopted. 
BZ integrations were carried out using Monkhorst-Pack 
sampling grids with resolution of 27rx0.04 for struc¬ 
ture optimizations. The atomic positions and lattice 
constants were optimized using the conjugate gradients 
(GG) scheme until the force components on each atom 
were less than 0.01 eV/A. An adequate fc-point sampling 
(30x30x 1) was employed for charge density distributions 
and band structure calculations. Phonon calculations by 
supercell approach as implemented in the PHONON code 
[29] was used to examine dynamical stability, and first- 
principles molecular dynamics simulations under con¬ 
stant temperature and volume (NVT) were performed 
to check thermal stability. 

During the structure searches, graphene always oc¬ 
cupies the lowest position in the enthalpy evolutions. 
We present a typical example (20 atoms/cell) in Fig. 
1(a). Besides graphene, phagraphene and <5-graphyne 
[10] with a 20-atom primitive cell were generated dur¬ 
ing this search. As we can see, the energies for graphene, 
phagraphene and (5-graphyne are -9.23, -9.03 and -8.49 
eV/atom, respectively. With different atom numbers in 
cell, many 2D carbon structures with hybrid sp- and sp^- 
carbon atoms were generated in our searches, but these 
structures usually have higher energies than sp^-carbon 
allotropes. For sp^-carbon structures, many of them are 
composed of penta-, hexa-, hepta-, and octa- (5-, 6-, 7-, 
and 8- ) carbon rings. Most of the planar carbon al¬ 
lotropes proposed in previous works can be reproduced 
in our systematic searches. Their energies and planar 
atomic densities are plotted in Fig. 1(b). It is noted 
that Dirac allotropes are mostly at much higher energies 
compared with graphene (denoted as 1), and structures 
with higher energies usually have lower planar atomic 
densities. The planar atomic density of phagraphene (de¬ 
noted as 27) is 0.37 atoms/A^, slightly smaller than 0.38 
atoms/A^ of graphene. Thanks to this dense packing of 
sp^-carbon atoms, it has lower energy than most of before 
predicted 2D carbon allotropes. 

Figure 2(a) shows its primitive cell and the periodic 
structure with notable space-inversion symmetry. The 
structure of phagraphene is composed of 5-6-7 carbon 
rings. Unlike P6m of graphene and i5-graphyne, the plane 
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FIG. 1. (color online) (a) Typical enthalpy evolution for 
a 20-atoni carbon system during an evolutionary structure 
search; (b) Energies and planar atomic densities for most of 
2D carbon allotropes re-obtained from our evolutionary struc¬ 
ture searches [30]. Allotropes with Dirac cones are labeled as 
blue triangles. 


group of phagraphene is Pmg. Being surrounded by 6- 
carbon rings, the pairs with 5- and 7-carbon rings invert 
to each other and arrange in a rectangular sublattices. 
Fig. 2(a) also presents the inequivalent carbon atoms 
(G1-G6) in its unit cell. It is noted that the length of Gl- 
C2 covalent bond is about 1.52 A, which is longer than 
that of G2-C3 (1.44 A), G3-C4 (1.41 A), G4-G5 (1.43 A), 
and C5-G6 (1.40 A) from our first principle calculations. 
Moreover, based on our band structure calculations as 
shown in Fig. 2(b) and 2(c), the valence and conduction 
energy bands meet at the Fermi level and show the fea¬ 
tures of a distorted Dirac cone in its first BZ [elaborated 
in Fig. 2(b)]. The Fermi velocity {vf) of phagraphene 
in both kx and ky directions were obtained by fitting 
these two bands at ki=Ki-\-q {i = x, y) to the expression 
of / = E{q)/h\q\, which can be evaluated via the slope 
of the bands. In the kx direction, dE/dkx =26.8 eVA 
(u/a;=6.48x 10®m/s); while in the ky direction, the slope 
of the bands equals -25.8 eVA (u/y=6.24x 10®m/s) and 
I4.2eVA(u/j^=3.43xl0’5 m/s). These Fermi velocities are 
comparable with those of other carbon Dirac allotropes 
predicted in previous works, as listed in Table I. 
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FIG. 2. (color online) (a) Structure of pliagraphene with notable space-inversion symmetry, C1-C6 are inequivalent carbon 
atoms in its unit cell; (b) The distorted Dirac cone formed by the valence and conduction bands in the vicinity of the Dirac 
point; (c) The comparison of band strnctures from density functional theory (DFT) and tight binding (TB) model. Inserted 
first Brillouin zone (BZ) with high-symmetry k points: F (0,0,0), X (0.5,0,0), Z (0.5,0.5,0), Y (0,0,5,0) and Dirac point: D (0, 
0.377,0); (d) Charge density distributions near the distorted Dirac cone, both Dirac bands (denoted as I and II) are from pz 
orbitals of sp^-carbon atoms. Fermi level has been set to zero. 


TABLE I. Calculated total energy (£'t)(eV/atom), planar 
carbon density (D)(atom/A^), plane group (G) and Fermi ve¬ 
locities (u/)(xl0®m/s) of the most of carbon allotropes with 
Dirac cones from GGA-PBE results. 


“ 

Et 

D 

G 


o-graphyne 

-8.30 

0.19 

P6m 

6.77 

/3-graphyne 

-8.38 

0.23 

P6m 

3.87/4.35/6.77 

5-graphyne 

-8.49 

0.26 

P6m 

6.96 

6,6,12-graphyne 

-8.51 

0.27 

Pmm f 

).56/6.04/6.29(Cone I) 





1.69/2.18(Cone II) 

phagraphene 

-9.03 

0.37 

Pmg 

3.43/6.24/6.48 

graphene 

-9.23 

0.38 

P6m 

8.22 


In order to explore the origin of the distorted Dirac 
cone, partial charge density distributions of the two Dirac 
bands in proximity of the Fermi level were first investi¬ 
gated, as shown in Fig. 2(d). It is clear that the oc¬ 
cupied states are inverted at the two sides of the Dirac 
point, which means the energy band inversion. Band-I 
and band-II, which are mainly contributed by the cou¬ 
pling of Pz orbitals along ai and 02 directions, cross each 
other and produce the Dirac point. In view of the cou¬ 
pling between pz orbitals leading to the formation of tt- 
conjugated framework, we adopted a simple tight-binding 


(TB) Hamiltonian of the 7r-electrons to describe the elec¬ 
tronic band structures in proximity of the Fermi level: 


H = -Y^ tijcf Cj + h.c. (1) 

b 


where tij is the hopping energy of an electron between 
the i-th and j-th atoms, c^^and Cj are the creation 
and annihilation operators, respectively. The distance- 
dependent hoping energy is determined using the for¬ 
mula: tij = to X exp[q X (1 — dij/do)], with to = 2.70 
eV, q = 2.20, and do = 1.5 A. Through diagonalization 
of a 20x20 matrix in the reciprocal space, the electronic 
band structures can be obtained, as indicated by the red 
dashed lines in Fig. 2(c). Obviously, this simple TB 
Hamiltonian in Eq. (1) can reproduce the band struc¬ 
ture of phagraphene in the vicinity of the Fermi level 
given by DFT calculations. 

As a matter of fact, most Dirac cones are indeed from 
the energy band crossing, which can be explained by a 
simple model: considering a 2-band system, the Hamil¬ 
tonian can be written as: 


H{k) = 


Hn{k)-E 

H2i{k) 


Hi2[k) 

H22{k)-E 


( 2 ) 
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FIG. 3. (color online) (a) Band structures of phagrahene 
along —r — Y — T line, with and without 5% strain along ai 
and 02 directions respectively. Fermi level has been set to 
zero; (b) Phonon dispersion of phagraphene. 


The appearance of Dirac cones corresponds to degenerate 
solutions of wave function of this Hamiltonian. In other 
word, the determinant of H(k) should be equal to zero 
and the following equations need to be satisfied: 

Hii{k) = H22ik);Hi2ik) = H2iik) = 0 (3) 

It is intrinsic that these three conditions in Eq. (3) must 
be simultaneously fulfilled to have a degeneracy, which 
is known as the von Neumann-Wigner theorem [31, 32]. 
In most cases, Hii{k) = H 22 {k) can be fulfilled by the 
space-time inversion symmetries, but the 2D system sym¬ 
metries cannot guarantee the existence of Dirac cones, 
because the number of variables {k^ and ky) is usually 
less than the number of equations to determine the Dirac 
points, which gives the reason of the rarity of 2D Dirac 
materials. As for phagraphene here, band-I and band-II 
in Fig. 2(d) or Fig. 3(a) can be simply considered as 
the two bands in above model, only they can appear in 
the vicinity of the Fermi level and invert to each other, 
the distorted Dirac cones are then produced at both sides 
of the space-symmetric and time-invariant fc-point Y. Al¬ 
though the determinant of its real Hamiltonian is hard to 
derive out analytical solutions, our numerical DFT and 
TB calculations all reveal that phagraphene happens to 
fulfill these three conditions simultaneously and presents 
distorted Dirac cones. 


Its direction-dependent character can be further inves¬ 
tigated by applying external strain, as shown in Fig. 3(a). 
The cone is stable under 5% strain along ai and a 2 di¬ 
rections, respectively. Due to the band symmetry, both 
position and anisotropy of the Dirac cone can be tuned 
by external strain. For 5% strain along ai direction, 
the Dirac point moves from (0, 0.377, 0) to (0, 0.417, 
0). The corresponding Fermi velocities change from 26.8 
eVA to 23.7 eVA in the k^ direction, from -14.2/25.8 eVA 
to -19.4/28.4 eVA in the ky direction. Things become 
different when a 5% strain is added along 02 direction. 
Dirac point shifts to (0, 0.313, 0), and the corresponding 
Fermi velocities become 28.1 eVA in the k^ direction and 
-25.4/6.7 eVA in the ky direction. 

In order to examine its dynamic stability, the phonon 
spectrum was calculated, as shown in Fig. 3(b), no imag¬ 
inary frequency modes were observed. To further exam¬ 
ine its thermal stability, a 3x3 supercell was built to 
perform ab initio molecular dynamics simulations. After 
being heated at room temperature (300 K) and 1000 K 
for 3 ps with a time step of 1 /s, no structural changes oc¬ 
cur, except some fluctuations in the 2D plane. All these 
indicate that phagraphene is a new low-energy carbon 
allotrope with distorted Dirac cones. Considering 5- and 
7-carbon ring dislocations and corresponding line defects 
have been observed in graphene [33], phagraphene may 
be realized experimentally in near future. 

In summary, systematic ab initio evolutionary struc¬ 
ture searches for 2D carbon networks identified a new 
low-energy 2D form of carbon, which is composed of 5- 
6-7 carbon rings and named as phagraphene. Thanks to 
its sp^-hybridization and density of atomic packing com¬ 
parable to graphene, this 2D carbon structure with Pmg 
plane group is lower in energy than most of the predicted 
2D carbon allotropes. Both DFT and TB model confirm 
its distorted Dirac cone in the first BZ. The direction- 
dependent Dirac cones are further proved to be robust 
against the external strain with tunable Fermi velocities. 
This exotic structure is not only promising for fully inves¬ 
tigating the massless Dirac fermions in 2D carbon elec¬ 
tronic systems, but also helpful for understanding the 
topological and correlated phases in the corresponding 
photonic artificial lattices. 
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